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ABSTRACT 

Context. Magnetic clouds (MCs) are formed by flux ropes (FRs) launched from the Sun as part of coronal mass ejections (CMEs). 
They carry away an important amount of magnetic flux and helicity. 

Aims. The main aim of this study is to quantify these quantities from in situ measurements of MCs at 1 AU. 

Methods. The fit of these data by a local FR model provides the axial magnetic field strength, the radius, the magnetic flux and the 
helicity per unit length along the FR axis. 

Results. We show that these quantities are statistically independent of the position along the FR axis. We then derive the generic 
shape and length of the FR axis from two sets of MCs. These results improve the estimation of magnetic helicity. Next, we evaluate 
the total magnetic flux and helicity crossing the sphere of radius of 1 AU, centered at the Sun, per year and during a solar cycle. We 
also include in the study two sets of small FRs which do not have all the typical characteristics of MCs. 

Conclusions. While small FRs are at least ten times more numerous than MCs, the magnetic flux and helicity are dominated by the 
contribution from the larger MCs. They carry in one year the magnetic flux of about 25 large active regions and the magnetic helicity 
of 200 of them. MCs carry away an amount of unsigned magnetic helicity comparable to the one estimated for the solar dynamo 
and the one measured in emerging active regions. Note: This article will be published in Solar Physics but is set in the A&A format 
because of LaTeX compilation problem with SP style. 
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1. Introduction 

Magnetic clouds (MCs) are a subset of Interplanetary Coro¬ 
nal Mass Ejections (ICMEs) characterised by an enhanced and 
smooth magnetic field strength, a large and coherent rotation 
of the magnetic field and a low proton temperature compared 
with the typical solar wind with the same velocity ( [Burlagaef a/.| 
1981| l. They are the continuation in the interplanetary medium 
of CMEs launched from the solar corona after an instability has 
occurred in the coronal magnetic field. Because of its observed 
properties, the large-scale magnetic configuration of MCs is fre¬ 
quently modeled by a magnetic flux rope (ER). 

Among others, two global quantities characterise a flux rope: 
its axial magnetic flux F and its magnetic helicity FI. This last 
quantity quantifies how all the elementary magnetic flux tubes 
are winded around each other in a defined volume. Magnetic 
helicity has several remarkable properties both from the theoret¬ 
ical and observational view points {e.g., see the reviews of |De-| 
|moulin|[2007t |Demoulin and Pariat||2009t [Pevtsov et a/.||2014| l. 

In particular, H is an ideal magnetohydrodynamic invariant that 
can be obtained from an invariant associated with electrons in a 
proton-electron multifluid description, in the limit of zero elec¬ 
tron inertia (see e.g., [Andres et q/.|| 2014} . In a closed system, 
magnetic helicity is almost conserved in resistive MHD on a 


largely transformed into other forms of energies. This theoret¬ 
ical prediction was tested positively with MHD simulations of 
coronal jets ( [Pariat et a(.|2015| l. 

The axial magnetic flux F and the magnetic helicity FI are 
conserved during the PR propagation unless the PR significantly 
reconnects with the surrounding solar-wind magnetic field. This 
conservation property was used to quantitatively link PRs ob 
served in situ to their solar sources (e.g., [Dasso eta(.|2005al[Lu 


|2005at 

oni et a/.|2005[|Mandrini et al. 2005 [ Qiu et a/.|2007[ Rodnguez 

et al. |20()8[[Hu et a/.|2014| l and torelate the in situ observations 

of two spacecraft, at 1 and 5.4 AU, of the same MC ( [Nakwacktj 
\et al.\20\\\ . Quantification of H and F also allow us to constrain 
models of coronal formation and ejection of flux ropes to the in¬ 
terplanetary medium, as well as of the dynamical evolution of 
MCs in the solar wind (for a review see [Dasso|2009[ ). 

The computations of F and even more of H are challenging 
because magnetic data are only available along the spacecraft 
trajectory, so along a ID cut of the PR, while these global quan¬ 
tities are 2D and 3D, i.e. they are surface and volume integrals, 
respectively. Then, their estimation relies on flux rope models 
with the free parameters of the model typically determined by a 
least square fit to the in situ data (e.g., [Al-Haddad et a/.|2013| 
references therein, and Section 2.3 1 . All the models provide an 


timescale lower than the global diffusion timescale (Matthaeus 
[and Goldstein|198^|Berger|1984[), while e.g. magnetic energy is 


estimation of the magnetic field within a cross section of the PR, 
so they provide F as well as H per unit length along the axis. 
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Then, H can be estimated with an ad hoc length of the FR, which 
is typically in the range [0.5,2.5] AU for a MC observed at 1 AU 
(see Section 2.3 i. This supposes a FR uniformly twisted along 
its length. 

Moreover, some theoretical models have been proposed to 
describe the global shape of MCs. They have been compared 
with in situ observations made by a single spacecraft {e.g., 
Marubashi and Lepping|[2()07[ [Hidalgo and Nieves-Chinchilla| 


2012|l. But the reconstruction of the 3D global MC shape from 


in situ measurements of a single event is not satisfactory because 
it is an ill-posed problem with no unique solution, and these 
models contain so many free parameters that generally several 
solutions compatible with the observations are found. In a re¬ 
cent work, [Janvier, Demoulin, and Dasso| ( |2013[ l proposed a new 
method to analyse the global shape of the main axis of MCs, 
from the statistical distribution of the orientation of a large sam¬ 
ple of events. This method was recently used to compare ob¬ 
servations with different models for the shape of the MC axis 
Panvier et a/.|2015] l. They found that an ellipsoidal shape is the 
model that best fits the data, and got an aspect ratio of ~ 1.2 for 
the ellipse. 

The main aim of this study is to further develop the compu¬ 
tation of H for MCs based on a statistical analysis of two sets 
of MCs. In Section]^ we first summarise the type of data used, 
then the equations needed to derive F and H, and finally we 
summarise our present knowledge of H estimations in MCs. In 
Section]^ we investigate how the main flux rope parameters are 
function of the curvilinear abscissa along the FR axis. In Sec¬ 
tion 1^ we propose a new method to estimate the length of MC 
axis; it is based on a statistical study of two MC sets. We next 
use in Section]^ the results of the previous sections to derive the 
amount of flux and helicity launched from the Sun per year and 
over a solar cycle by MCs/CMEs. We compare these results with 
the contribution provided by the much more numerous small FRs 
detected in the solar wind at smaller scales than MCs, as well as 
with other solar estimations of magnetic helicities {e.g. dynamo, 
emerging active regions, solar wind). Finally, in Section]^ we 
summarise our results and outline future studies needed to im¬ 
prove the global helicity budget. 


2. Observations and Models 


2.1. Data Sets 

In order to perform a statistical study of the magnetic flux and 
helicity of MCs, we select the two largest lists of analysed MCs 


presently available. Lynch et al. ( |2005| l studied 132 MCs ob¬ 
served nearby Earth by Wind and ACE spacecraft during the 
period 1995-2003. [Lepping and Wu| ( |2010| l studied 98 MCs 
observed by Wind spacecraft. This list was extended to the 
time period of Eebruary 1995 to December 2009 (Table 2 at 
http;//wind.nasa.gov/mfi/mag_cloud_Sl.html Removing a few 
MCs that were badly observed (crossing too close from the 
boundary), it remains 107 MCs (see [Janvier, Demoulin, and] 
|Dasso|[2013| for more information). Below, we refer to MCs 
of both lists as MCLy and MClb sets, respectively. 

The local magnetic configuration of the studied MCs was de¬ 
duced in both studies by following the fitting procedure of jLep-j 
[ping, Burlaga, and Jones| ( [T990| l, i.e. with a least square fit of the 
magnetic held data along the spacecraft trajectory with a linear 
force-free magnetic held having a circular section and a straight 
axis ( Lundquist|| 1950 1. The linear force-free held corresponds 
to the relaxed state with minimum energy for a given helicity 
content and axial held distribution. In the ER coordinates, with 
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Fig. 1. Schema showing the definition and the large-scale meaning of 
the location angle A for a FR launched from the Sun. The FR structure is 
outlined by one twisted field line (black) with dotted style for the behind 
side. The schema shows the plane of the FR axis which is inclined by 
an angle i on the ecliptic. A is defined by the angle between the radial 
direction [Up] and the normal to the axis [n]. The cylindrical coordinates 
of a point along the axis are (p, (p). The full range of is 2 (^Pmax. The 
signed curvilinear coordinate [j] is defined along the FR axis with its 
origin set at the apex. 


z along the FR axis, the magnetic field Bl of the Lundquist’s 
model writes 


Bl = Bo[7i(a'r)ea -h JQ{ar)tj] , 


( 1 ) 


where Jq and Ji are the ordinary Bessel functions of order 0 
and 1, and Ca and Cz are the azimuthal and axial unit vectors in 
cylindrical coordinates. Bq is the magnetic field strength on the 
axis and a is the linear force-free constant. The authors selected 
the boundaries of the MC such that the magnetic field becomes 
purely azimuthal there, i.e. they selected |a| = cjR, where the 
constant c is the first zero of Jq (c a; 2.4) and R is the FR ra¬ 
dius. The handedness of the FR is defined by a another param¬ 
eter (which values are +1). It is equivalent to define a signed a 
parameter. 

The orientation of the axis is defined by its longitude [0] 
and latitude [6] in the Geocentric Solar Ecliptic (GSE) system 
of reference. Another parameter is the closest approach distance 
[To] between the spacecraft trajectory and the MC axis. It is 
frequently normalised to the FR radius, and called the impact 
parameter [p - Tq/R]. Taking into account the observed mean 
velocity inferred from the in situ plasma measurements, the least 
square fit of the above model to the magnetic data determines 
the six parameters: Bq, R, (f), 6, p and sign(a). The fit is done in 
two steps. In the first step, each magnetic field vector is divided 
by its norm (to avoid a bias due to a typical asymmetry of field 
strength between the front and rear of MCs). In the second step, 
a fit to the data is realised to determine Bq (keeping the other 
parameters fixed). The quality of the MC fit is me asured by the 
square root of the chi-squared Xiu - V?/^ for Lynch et al. 
( |2005 I and by the square root of the reduced chi-squ ared defined 
as xr - -\jx^l0^d ~ ”) for Lepping and Wu ( 2010| where n - 5 
is the number of parameters of the fit and Ad is the number of 
data. Both Xdir and are computed during the first step of the 
fit: they are dimensionless quantities, and measure how well the 
model fits the direction of the observed magnetic field. 
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We also analyse two other lists of interplanetary FRs. They 
can be found in the papers of [Feng, Wu, and Chao ( 2007| l and 
Feng et al. (2008|l. The detected FRs were also fitted with the 


Lundquist field with a similar procedure as [Lepping, Burlagai] 
and Jones ( 1990|l, and the derived list of events contains mostly 


f B^(r') 2n / dr' = BqR^ , 

Jo c 


( 2 ) 


H^L f 

Jo 


2AJr) B^(r) 2nrdr - 


27iJ]{c) 


bIfJl . 


(3) 


For a FR configuration, magnetic helicity is directly related 
to the mean number of turns per unit length {nj of the magnetic 
field lines along the axis (Section]^. More precisely: 


the same parameters (see the above papers for the small differ¬ 
ences). The list of |Feng, Wu, and Chao ( |2007[) has 144 FR s with 
both MCs and small FRs, while the list of Feng et aT] p008| l 
contains 125 small FRs. 

[Janvier, Demoulin, and Dasso| ( |2014al l have analysed the dis¬ 
tributions in MC radius of the four data sets and found that small 
FRs and MCs have different distributions: a power-law for small 
FRs and a Gaussian-like distribution for MCs. They concluded 
that small FRs have different solar origins compared with MCs. 
We refer to this study for further information on the data sets. 

[Janvier, Demoulin, and Dasso| ( |2013] l have introduced a new 
spherical coordinate system for the FR axis since the longitude 
(j) has a large error for FR oriented close to the north-south direc¬ 
tion. They set the polar axis of the new spherical coordinate sys¬ 
tem along Xgs£, then they defined the inclination on the ecliptic 
[i] and the location [T] angles. A is defined by the angle between 
the radial from the Sun and the normal to the axis (Figure [T). For 
a FR with a known axis shape, of the type shown in Figure [11 the 
angle A has a monotonous variation along the axis. For these 
cases, A labels the location along the axis where the spacecraft 
intercepts the FR. The sign convention of A is such that, for a 
FR close to the ecliptic (i small) the eastern (western) leg, corre¬ 
sponds to T < 0 (T > 0). This sign convention is extended to all 
i values by continuity. 

2.2. Theoretical Estimations of Global Quantities 

Two main global quantities of a FR are its axial magnetic flux [F] 
and its magnetic helicity [//]. These global magnetic quantities 
can only be estimated from the fit of in situ data by a FR model. 
H is first estimated per unit of length along the FR axis, and it is 
typically given for a fixed length, which values are discussed in 
Section|2^ 

Below we write the expressions of F and FI for a magnetic 
field having a cylindrical symmetry, a local approximation for 
the FR of MCs. Then, B(r) = BJr) + B^(r) Cz, where B^, B^ 
are the azimuthal and axial components depending only on the 
radial coordinate r. Next, we write the specific results for the 
Lundquist model (Equation ([1])). 

The axial flux, integrated from the axis to the FR radius R 
and assuming a cylindrical symmetry, is 


where the constant c is the first zero of the Bessel function Jo(r). 

The relative self magnetic helicity of a flux rope is the sum 
of its twist and writhe helicities ( |Berger and Prior||2006 1. For 
MCs it is mostly limited to the twisted helicity since the FR axis 
is thought to have a low writhe as shown for a few MCs ob- 


H ^ ill F^ L, 


(4) 


where is a flux weighted mean of the number of turns per unit 
length along the axis. 


Ht 




(5) 


where n{F') is the local number of turn, which can be expressed 
in term of the cumulated flux F'{r) from the origin of the flux 
rope up to a radius r, as shown in Section and where the 
right-hand side expression is for the Lundquist model. 


2.3. Helicity Estimation of Magnetic Clouds 


Since magnetic helicity is intrinsically a 3D quantity (i.e., com¬ 
puted with a volume integral) while observations are limited to 
the magnetic field measured along the spacecraft trajectory, the 
estimation of MC helicity involves hypotheses and models. 

Within the cylindrical hypothesis, different MC models have 
been proposed. For example, a uniformly twisted field (Dasso 
et al 2003 , a non-force-free field with constant current HE] 
dalgo et al. |2000|l or with an azimuthal component of the cur¬ 


rent depending linearly on the radius ( Cid et aT\2002\ . The fit 
of these models to data introduces a variation of the deduced 
helicity up to 30%, which still remains small compared to the 


variation of helicity computed between different MCs (Gulisano 
ef Q/.|20d5ll. Extensions to elliptical cross-section (e.g., [Vandas 


and Romashets|p003|l increase the helicity approximately pro¬ 


portionally with the aspect ratio of the cross section (|Demoulm 


served by several spacecraft {e.g., Burlaga, Lepping, and Jones 
|1990[|Ruffenach ef af.|2012] l. In terms of order of magnitude, the 
writhe contribution is of the order of 0.1 equivalent turn or be¬ 
low, while the twist is important, of the order of 10 turns. Then, 
one can consider that the MC helicity is mostly due to the twist. 
The helicity [//] of a straight flux rope of length L is ( |Berger| 
|2003[[Dasso et a/.|2003| l: 


land Dass'o||2009b|l. With typical values o f 2 to 3 of the aspect 
ratio ( [Demoulin, Dasso, and Janvier|2013| l, this increases signif¬ 
icantly the estimation of the helicity values. Einally, non force- 
free models {e.g.. Mulligan a nd Russell||200l] |Hidalgo]|201 1[ 
[Isavnin, Kilpua, and Koskinen|201 l| l have also been developed. 
It would be worth to develop both their helicity estimations and 
their applications to a larger number of MCs using these models, 
but it is out of the scope of the current paper. 

The above models have a straight axis configuration. Then, 
the derived helicity is only a local estimation per unit length 
along the axis. An appealing approach is to extend the above 
models to toroidal geometry in order to include the curvature of 
the ER axis {e.g., |Marubashi and Lepping|2007{|Romashets and| 
[Vandas [2009| l. Due to a larger number of free parameters, it is 
not yet clear if they can all be constrained by the data of a single 
spacecraft. Two well separated spacecraft provide more con¬ 
straints to the toroidal model ( [Nakagawa and Matsuoka|[2010| l. 
However, the number of MCs observed is very limited in this 
configuration as it requires a ER oriented close to the ecliptic 
plane, where spacecraft are typically located (see the review of 
Kilpua et al. 201 1[ ). These models also assume an invariance 
along the curved axis. In fact, it is not known how the twist 
is distributed along the MC axis. We estimate this dependence 
from a statistical study in Section]^ 

Most MCs are faster than the local solar wind, at least close 
to the Sun. It results in the formation of a sheath before the MC 
where plasma and magnetic field accumulate. When magnetic 
fields of different orientation are pushed together, it generally 
implies magnetic reconnection. This phenomenon can also oc¬ 
cur at the rear of the MC, for example, when a faster MC or a 
fast solar wind stream takes over the propagating MC. This leads 
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to a FR progressively peeling off and only the central region re¬ 
mains as a coherent FR when observed (|Dasso et al.|2006 1 . This 


is confirmed by the presence of magnetic discontinuities (Dasso 


\et a/.|[2007 [Nakwacki et al.\20\\) and by in situ reconnection 
signatures ( Ruffenach et a/.|2012| 2015] l. The amount of recon¬ 
nected flux is case-dependent and significant: it was found that 
about 40 % of the total azimuthal magnetic flux on average is 
lost from this erosion process. This process of erosion is typi¬ 
cally not taken into account in the definition of MC boundaries, 
in particular for the lists of MCs of [Lynch et af.| ( [2005| l and |Lep-| 
ping and Wu (2010|l. As such, the helicities estimated in the 


present paper are expected to be intermediate between those of 
the FR remaining at 1 AU and of the FR before erosion. 

The reconnection of the FR with open solar wind field can 
also occur in one leg without direct consequences for the in situ 
magnetic measurements (when the magnetic perturbation has no 
time to travel to the crossing location). However, the tails of the 
electron distributions provide clues on the large scale connec¬ 
tivities since the faster electrons are fewly interacting with the 
plasma. The presence of bi-directional, or counter-streaming, 
electron heat fluxes in a MC is generally interpreted as a con¬ 
nection to the Sun at both field-line ends ( [Richardson, Farrugia, 
and Burlaga|199l[|Shodhan et a/.|2000|l. The counter-streaming 


electrons are typically present in fragmented and partial (from 0 
to 100%) portions of MCs observed at 1 AU ( [Shodhan et ai 
2000 [l. The counter-streaming electrons in MCs observed at 


about 5 AU are also very case-dependent, and they are present 
on average 55% of the time, a result comparable to observations 
at 1 AU, so that the amount of disconnection from the Sun does 
not increase with distance ( [Crooker et a(.|20()4| . 

The length of field lines can be inferred in exceptional cases, 
when high-energy electrons, accelerated close to the Sun, are in¬ 
jected in them and detected in situ. A first method is to derive the 
path length from the velocities and the different arrival times of 
electrons of various energies ( [Kahler and Ragot][2006[ [Mass on | 
\et a/.|20f^ . The second method is based directly on the travel 
time, so it requires an estimation of the solar release time (with 
the onset of type III radio bursts) and the in situ detection of the 
same electron beam. The unique MC analysed by Larson et al. 
\\991\ was recently extended up to a list of 18 MCs ([Kahler, 


Krucker, and Szabo|201 l[[Hu et q/.|2014[[Hu, Qiu, and Krucker 


2015[ l. The last studies show an estimation of the length across 

the flux ropes consistent with what is expected for a flux rope 
with a more uniformly distributed twist across the cross-section 
than what is predicted by the Lundquist’s model. All these stud¬ 
ies give a large range of lengths, from below 1 to 4 AU. Since 
only a few field lines can be probed in one FR leg at most, esti¬ 
mating the flux rope axis length is therefore rather limited. These 
results would need to be interpreted with an estimation of the lo¬ 
cation of the spacecraft crossing (so of the values of the location 
angle A and impact parameter p), and with the knowledge of 
which FR leg the energetic electrons were traveling in (i.e., the 
sign of T). Then, the estimation of the whole effective FR length 
is not straightforward from these results. 

In practice, for the magnetic helicity estimation of MCs, a 
length is typically assumed. Different values have been used, 
ranging from 0.5 to 2.5 AU ( DeVore[[2000tm'nch et a/.[[2005| l. 
For some specific cases, physical arguments have been invoked 
to justify the selected length such as the initiation of the so¬ 
lar ejection by the kink instability (Nindos, Zhang, and Zhang 
2003 , or the disappearance of the solar source region ( [Mandrim 
et al. 2005|l, or the agreement between the azimuthal flux esti¬ 


mated in the MC and the flux swept by the flare ribbons in the 
source region ( [Du, Wang, and Hu|2007[ [Hu et al. ]20T^ . 
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In conclusion, the axis length of the FR is still a major source 
of uncertainty to estimate MC helicity. Furthermore, it is not 
known how helicity is distributed along the MC axis. The aim of 
the next sections is to improve our knowledge on these aspects. 


3. Dependence Along the Flux Rope Axis 

In this section we test whether the MC properties are variable 
along the MC axis, so whether there is a statistical dependence 
with the location angle A since this angle is also a coordinate 
along the axis with |/l| increasing away from the apex (Figure[^. 
We characterize the correlations by two coefficients: the Pearson 
[cp] and the Spearman rank [cs] correlation coefficients. We re¬ 
port in different figures the fit of the data by a linear function to 
show global tendencies, as well as the mean value of the studied 
property [/r] and the standard deviation of the fit residuals [cr] 
(computed with respect to the fitted straight line). 

3.1. Influence of the Spacecraft Trajectory Location 

The parameters Xdn and xr are both testing how well the mag¬ 
netic field direction of the model fits the in situ data (Sec- 
tion[2.1|i. For the MCpy set, there is a weak tendency (cp = cs = 
-O.lh) of a lower ;t'dii as |/1| increases for both MC legs (Fig¬ 
ure [^). This tendency becomes even weaker when the MCs, 
crossed near their outer boundaries, are removed from the sam¬ 
ple. For example, cp = cs = -0.07 with the selection \p\ < 0.7 
where p is the impact parameter. For the MCpe set there is no 
significant correlation;^fR(/l) for both MC legs (Figure]^). This 
result is robust as it is also valid for sub-groups of MCs {e.g. 
Cp = 0.04, Cs = 0.06 with the selection \p\ < 0.7) and there is 
no significant differences between both legs. We conclude that 
the quality of the Lundquist fit to the data is independent of the 
spacecraft crossing location along the flux rope. 

The impact parameter p is spread in the interval [0,1] as 
expected with random distance encounters. Still, low p val¬ 
ues are significantly more numerous (Figure [^. This is a con¬ 
sequence of the oblateness of the flux-rope cross section s 
[moulin, Dasso, and Janvier[2013] l. The correlations of the impact 
parameter p with A are positive for MCpy and MCpe sets (Fig- 
ure[^,b) and all correlation coefficients are small with the selec¬ 
tion |/1| < 50° (|cp|, |cs| < 0.04 for MCpy set and |cp|, |cs| = 0.01 
for MCpe set. Figure [^,d). We interpret this change of p in 
the flux rope legs with an observational bias, as follows. As the 
spacecraft trajectory is crossing the flux rope less perpendicular 
to its axis (larger |/1|), the spacecraft trajectory explores a longer 
part along the flux rope. There, the bending of the MC axis af¬ 
fects the measurements of the magnetic held. It implies that the 
hypothesis of a local straight flux rope, used in the Lundquist 
model, is less valid as |/1| increases ( [Owens et a/.[|2012l l. In¬ 
deed, for moderate |/1| values, the deviation between a curved 
and straight flux rope is small but it becomes strong for large 
|/1| values (see their Figures 3-5). This deviation is interpreted 
by the Lundquist’s fit as a larger impact parameter so that p is 
positively correlated with |/1| in Figure[^,b. 

3.2. Variation of the Physicai Parameters Aiong the Axis 

For both MCLy and MClb sets there is no significant correlations 
between the axial field strength Bq with A (Figure |^,b). This 
result is robust as it stays for sub-groups of MCs. For example, 
Cp = Cs = -0.1 for MCpy set and cp = -0.08, cs = -0.06 
for MCpe set with the selection \p\ < 0.7. It also implies that 
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Mq, -0.16, c,=-0.16 (a) 

Xd,T ^ = 0.13,0 = 0.07 



Mq, Cp= -0.05, c,=-0.10 (b) 

Xk pi = 0.14, 0 = 0.05 



Fig. 2 . Dependence of the fit quality measured hy xhy oi Xr (non dimensionalised) versus the absolute value of the location angle [d]. (a) and 
(b) are showing MCs analysed by Lynch et al. (2005) and Lepping and Wu (2010), respectively. The straight lines are linear fits to the data points 
showing the global tendency, d > 0 and d < 0 are respectively shown in red and blue, and the abscissa, |d|, allows the comparison of the two 
FR sides (Figuref^. The results with the full MC sets are shown in black, c,, and c, are respectively the Pearson and Spearman rank correlation 
coefficients, n is the mean value of the ordinate and cr is the standard deviation of the fit residuals. 



Fig. 3. Dependence of the impact parameter [p] versus the absolute value of the location angle [d] for the two sets of MCs: (a,c) MCcy and (b,d) 
MClc sets. The top row is for all MCs while the bottom row is for a reduced interval of |d| (< 50°). The straight lines are linear fits to the data 
points showing the global tendency, d > 0 and d < 0 are respectively shown in red and blue, and the abscissa, |d|, allows the comparison of the two 
FR sides (Figure [^. The results with the full MC sets are shown in black, q and are respectively the Pearson and Spearman rank correlation 
coefficients, p is the mean value of the ordinate and cr is the standard deviation of the fit residuals. 


Mq, Cp=-0.13, c,= 0.01 (a) 

^^0[nT] p=15nT, o = 7nT 



0 20 40 60 80 


Mq^ Cp= 0.03, c,=-0.04 /|a\ 

Bo[nT] p=17nT, a = 9nT 



Fig. 4 . Dependence of the axial field strength [Sq] versus |d|. (a) and (b) are showing MCs analysed by Lynch et al. (2005) and Lepping and Wu 
(2010), respectively. The straight lines are linear fits to the data points showing the global tendency, d > 0 and d < 0 are respectively shown in red 
and blue, and the abscissa, |d|, allows the comparison of the two FR sides (Figure[TJ. The results with the full MC sets are shown in black, q and 
Cj are respectively the Pearson and Spearman rank correlation coefficients, p is the mean value of the ordinate and cr is the standard deviation of 
the fit residuals. 
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Fig. 5. Dependence of the flux-rope radius [/?] versus the absolute value of the location angle [d] for the two sets of MCs: (a,c) MCcy and (b,d) 
MClc sets. The top row is for all MCs while the bottom row is for a reduced interval of |/l| (< 50°). The straight lines are linear fits to the data 
points showing the global tendency, d > 0 and d < 0 are respectively shown in red and blue, and the abscissa, |d|, allows the comparison of the two 
FR sides (Figuref^. The results with the full MC sets are shown in black, c,, and c, are respectively the Pearson and Spearman rank correlation 
coefficients, n is the mean value of the ordinate and cr is the standard deviation of the fit residuals. 





Fig. 6. Dependence of the number of turns per unit length [n,] versus the absolute value of the location angle [d] for the two sets of MCs: (a,c) 
MClj and (b,d) MClc sets. The top row is for all MCs while the bottom row is for a reduced interval of |d| (< 50°). The straight lines are linear fits 
to the data points showing the global tendency, d > 0 and d < 0 are respectively shown in red and blue, and the abscissa, |d|, allows the comparison 
of the two FR sides (Figure [Tl. The results with the full MC sets are shown in black, c,, and c., are respectively the Pearson and Spearman rank 
correlation coefficients, /j is me mean value of the ordinate and tr is the standard deviation of the fit residuals. 


when observed at 1 AU the axial field strength has no signifi¬ 
cant dependence along the MC axis. This result seems a priori 
contradictory to the standard picture of a MC {e.g. Figure 1 in 
Richardson and Cane|2010 1 , for which Bq would be stronger at 
the leg than at the apex of the MC. The result can then be un¬ 


derstood as follows. MCs strongly expand as they move away 
from the Sun as a consequence of the approximative balance of 
total (magnetic and plasma) pressure between the MC and the 
surrounding solar wind ( |Demoulin and Dasso|2009a i. This im¬ 
plies that Bo is mainly a function of the solar distance. Then, 
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we interpret the above uniform distribution of Bq along the MC 
axis as a consequence of an approximative pressure balance at 
a fixed observation distance (at 1 AU). Finally, the MCly and 
MClb sets have comparable Bq distributions, with a mean value 
and dispersion of 16 ± 9 nT. 


In contrast to Bq, R is statistically a decreasing function of 
|/i| (Figure l^,b) with a stronger anti-correlation for the MCly 
than for theMCLe set. This anti-correlation is due to the ab¬ 
sence of large R values for large |/1| as can be seen by the ab¬ 
sence of blue/red dots in Figure [^,b. Indeed, with the selection 
|/1| < 50°, the correlations are much weaker for both sets of MCs 
(Figure ^,d). The correlations are even weaker if the stronger 
criteria /I| < 40° is applied (cp = -0.001, cs = -0.04 for MCly 
set and cp = -0.07, cs = -0.06 forMCpe set). A selection on \p\ 
has a lower effect on the correlations. We interpret these lower 
values of R in the flux rope legs with an observational bias, as 
above for the impact parameter p. Next, the difference between 
both legs in Figure]^ is small and not confirmed by Figure]^. 
Finally, the MCpy and MCpe sets have comparable R distribu¬ 
tions with a mean value and dispersion of 0.12 + 0.05 AU. 


We have performed a comparable analysis with the mean 
number of turns per unit length [tit]. With a Lundquist model, 
with an axial field vanishing at the flux-rope boundary, is di¬ 
rectly related to R (Equation (^). Still, as nt(A) tells how the 
twist is distributed along the flux ropes axis, we also show the 
results with n^. With both the MCpy and MCpe sets, has a 
strong positive correlation with |/l| (Figure |^,b), which is sur¬ 
prising as the MC legs would be more twisted than the apex. In 
fact, with the selection |/l| < 50°, the correlations are weak for 
both sets of MCs (Figure |^,d). They are even weaker with the 
more stringent condition |/l| < 40° (cp = 0.09, cs = 0.04 for 
MCly ‘'P “ 0.01, Cs = 0.06 for MCpe set). We conclude 

that the flux ropes are uniformly twisted along their axis, at least 
in the range |/1| < 50° around the apex within the limits of the 
variations between MCs: =4 + 2 AU“'. The mean num¬ 

ber of turns found here is in agreement with previous studies. 
For example, [Farrugia et a(.| ( [T999[ ) studied a small and hot flux 
rope assuming a constant twist model (i.e., the Gold and Hoyle 
model) and found a number of turns of ~ 7 AU '. On the other 
hand, [Mdstl et a/.| ( |2009| l studied one MC from in situ observa¬ 
tions made with two spacecraft (STEREO and Wind) crossing 
different parts of the cloud, and modeling the magnetic topology 
with a Grad-Shafranov equilibrium. They found a small varia¬ 
tion of the number of turns across the flux rope, with a mean 
value ~ 2 AU”'. 

The above absence of significant correlation for Bq, R and 
Ht with A, at least for |/l| < 50°, implies also that the global 
quantities F and FI, Equations (|2][^, are also almost indepen¬ 
dent of 4. The correlations of F and H with A can also be done 
directly. However, the main limitation of this approach is the 
much broader range of variation within a MC set since non linear 
dependencies on Bq and R are present in F and H (Equations Q 
[^) and there is also a positive correlation between Bq and R 
(cp = 0.39, Cs = 0.35 for MCpy set and cp = 0.31, cs = 0.23 
for MClc set). This implies a much larger dispersion of these 
global quantities so that a correlation study is less pertinent {e.g. 
it is more affected by outsiders). This is especially true for H, 
which has the strongest non-linearities, while, in contrast, the 
mean number of turns [tii] has a relatively limited range of varia¬ 
tion within MCs, so we can better test its correlations (Eigure|^. 


4. An Estimation of the Flux Rope Axis Length 

Since we found no significant dependences along the ER axis 
of Bq, R and in the previous section, we simply need the ER 
length to estimate the total magnetic helicity of MCs. At 1 AU, 
this length was typically taken in the range [0.5,2.5] AU in previ¬ 
ous studies (Section [23 ] i. In this section, the length is estimated 
from the information derived statistically on large samples of 
MCs. 


4.1. Method to Derive a Mean Axis Shape 

With a set of MCs, one can define an observed probability dis¬ 
tribution for each parameter of the fitted model. Janvier, De- 
|moulin, and Dasso] ( |2013 1 have developed a method to deduce 
a generic MC axis shape from the observed probability !Pobs(4). 
The main idea is that MCs are crossed at various locations, i.e. 
at different A values along their axis. Then, the observed prob¬ 
ability !Pobs(d) is a consequence of the axis shape, with more 
detections expected as the local orientation of MC axis is further 
away from the radial direction from the Sun (Eigure[T]). The sta¬ 
tistical analysis supposes that all MCs have a comparable axis 
shape with only a scaling factor in the angular extension (2(f„iax, 
Eigure[^. Indeed, the probability !Pobs(^) was shown to be nearly 
independent of the MC parameters such as field strength, radius 
and inclination on the ecliptic when the MCs were analysed in 
sub-groups ( | Janvier, Demoulin, and Dasso||2013) . Eurther they 
showed with a synthetic MC axis model that the angular exten¬ 
sion 2^inax has almost no effect on !Pobs(4). This justifies the 
analysis of all the MC set together, and the derivation of a mean 
axis shape from !Pobs(d) with (.Pmax as the only free parameter. 


Below, we first shortly summarise the analysis of Janvier, 
|Demoulin, and Dasso] ( |2013| l before extending it to derive the 
curvilinear abscissa along the axis, and then its length. The 
MC axis is supposed to be inside a plane and it is described 
with cylindrical coordinates \p, (f\ (Eigure [^l. The probability 
of crossing a MC can be expressed either in function of ip as 
V^piipldip, or in function of A as !Pobs(4) dA and these two proba¬ 
bilities are equal. Since CMEs are launched from a broad range 
of solar latitude and any longitude over the time scale of the anal¬ 
ysed MC set (almost a solar cycle), the MCs are expected to be 
crossed with a uniform distribution in tp, so = 1 /(2^n,ax) with 
the normalisation of the probability to unity. At the difference 
of [Janvier, Demoulin, and Dasso] ( [2013| ), here we do not sym¬ 
metrise P^((p) and Pohsid), so we keep separated positive and 
negative values of ip and A. The above equality of probabilities 
implies 


d(p — 2 (^rnax P obs(4) dA . 

Its integration defines ^ as a function of A as 


<^max) — 7 (pjj 


■f" 


Pob.(d') dA'. 


( 6 ) 


(7) 


Next, we relate p to 4 by expressing 4 as the angle between the 
radial direction [Up] and the normal to the axis [n] (Eigure 
which writes as 

d Inp = - tan(4) d^. (8) 

Using Equation (|^, the integration of Equation ([^ implies 

lnp(4, (^max) — “7 ipjjiax I tan(4 ) P obs(^ ) ti4 + Inpjnax ■ (9) 

Jo 
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Equations (0 and (|^ define a generic flux rope shape as a para¬ 
metric curve (p{A),ip{X)) in cylindrical coordinates, from the 
probability distribution !Pobs(d) derived from the analysis of a 
MC set. 

We extend the previous analysis by defining the curvilinear 
elementary length [ds] along the axis as 

d5 = 7(dp)2 + (p d^)2 = ^1 + j P , (10) 

after introducing Equation Proceeding as above for the 
derivation of (p{X) and p(A), the curvilinear abscissa with origin 
at the apex is 

rV(^')nbs(^') 

^max) “ 2 ^max I /i/\ ‘ (H) 

Jo cos(/l') 

At the limit |/1'| ^ 90°, cos(/l') ^ 0 at the denominator. 
However, in the observations !Pobs(d') strongly decreases with 
|/i'| and vanishes above |/i'| > 80° (Eiguresj^ and|^) so that the 
integral is not singular, but rather s{A) is typically flat for large 
|/i| values (Eigures|^ and[^). 


4.2. Mean Axis Length 

The application of Equations (|^ and Q to the MClc set is shown 
in Eigure|7]r for three values of The deduced axis shape is 
only weakly asymmetrical between the two sides, then compara¬ 
ble to the schema drawn in Eigure 2 of |Burlaga et al. ( 1998|l. This 
result is also comparable to Eigure 2 of Zurbuchen and Richard- 1 
son] (2006 1 for the front part while the legs are bended with a 
Parker-like spiral. The spiral shape is due to the rotation of the 
Sun carrying the anchored field lines. 

The curvilinear abscissa is evolving more around the apex 
(A = 0) than at larger |/l| values (Eiguresj^ and[^). There is 
also a weak asymmetry between the legs as shown by compar¬ 
ing the continuous curves with the dashed curved showing s(/l) 
computed with a symmetric probability (imposing !Pobs(d) = 
Pobs(-'i))- As expected, the angular extension has a sig¬ 
nificant effect on the curvilinear abscissa (Eigure |^). 

The results with the MCby set are close to the ones for MClb 
set with the minor difference of s(/l) increasing more sharply 
close to T = -1-80° for MCby set (Eigure |^). This is a con¬ 
sequence of the local maximum in the positive tail of !Pobs(d) 
(red curve in Eigure [^). It has also implications for the axis 
shape with a positive A leg extending more towards the Sun for 
MCby than MClb set (Eigure|^). The other local peaks or dips in 
^ob.s(T) (Eigure 1^) have only a weak effect on s(/l) (Eigure [^). 
This is a consequence of the integration averaging effect. 

The quantities Bq, R and tit are statistically independent of A 
in the range [-50°,50°] (Section [3. 2| l. Then, a minimal length 
to compute H is within this A interval (unless the flux rope re¬ 
connects with the solar wind magnetic field). Another estimate 
of the length is to extend this interval to the full range of d; 
[-90°, 90°]. More generally, the length can be estimated for the 
range [-dsup, dsup] as 


Tint(dsup5 ‘fmax max )- s(- 

■'^sup^ ^max) ■ 


( 12 ) 


Eigure shows the evolution of Lint(dsup, ^f>max) in function of 
ijOmax- It IS nearly a linear function of i^Cmax, because the curvilin¬ 
ear abscissa s(/lsiip, ^max) is proportional to ^max (Equation (111). 


However, the linearity is only approximative because an extra 
dependence on (f„iax is present in p{A') (Equation (|^). This de¬ 
pendence is weaker since Inp, and not p, is a linear function 
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of ifimax- Next, there is only a slight increase of Lim between 
Tsup = 50° and dsup = 90° (thin and thick lines, respectively, 
in Figure ^). The lengths computed with the two different sets 
MCby and MClc are comparable (red and black continuous lines, 
respectively). The slightly larger Lint for djup = 90° and MCpy 
set is a consequence of the larger extension of the computed 
axis towards the Sun (Figure |^). This effect is reversed for 
Tsup = 50°. Finally, for around 30°, the typical CME ex¬ 
tension observed with imagers (see below). Lint could simply be 
approximated by the linear function 0.2 -1-3.2 ^max/90. 

If the flux rope is still attached to the Sun by both legs, 
a lower estimate of the total length is given by adding radial 
straight lines linking the photosphere to the ends of the axis 
shape found above (Figures]^ and [^. Since the Parker spiral is 
close to the radial direction close to the Sun, this straight line ap¬ 
proximation provides only a slight underestimation of the length. 
This total length writes 


Ttotal(Tsup, l^max) — Tint(dsup, *^max) 

+ p(Tsiip, ^Cmax) + P(~d sup? ^max )-2f?o-(13) 

where Rq is the Sun radius. Ltotai is even closer to a linear func¬ 
tion than Lint (Figure]^), because the contribution of the straight 
leg parts (p(+/lsup, ipmax) - Rq) nearly compensates the contribu¬ 
tion of Lint with increasing ^max values. The same is true for the 
dependence with Tsup: Ltotai(Tsup, ^max) curves superpose each 
other very well in Figure for Asup = 50° and 90°. Finally, the 
linear function 2 -i- 1.8 ^max/90 approximates Ltotai very well. 

An estimation of (f„iax with in situ measurements is gener¬ 
ally not possible because only a few MCs are crossed by several 
spacecraft (e.g., Burlaga, Lepping, and Jones|1990[[Ruflenach 


et al. 20121. However, its estimation could be given from CME- 


imaging, which records a higher number of CMEs. Observations 
of CMEs situated close to the Sun and at the limb minimise the 
projection effects, although the tilt of the flux rope axis cannot be 
inferred. Since the orientation of the flux rope is not determined, 
this supposes a comparable angular extension of CMEs along 
and across the flux rope. The typical nearly circular observed 
shape of CMEs directed toward the observer (full halo CMEs) 
justifies this hypothesis. Wang et al. ( 201 1[ ) derived from the 
Large Angle and Spectrometric Coronagraph (LASCO) a mean 
(Pmax of 30° for limb CMEs and for 65% of the CMEs, its val¬ 
ues lie within the interval [15°,45°]. Since these values are de¬ 
duced from coronagraphs, which image the densest parts of the 
CMEs, namely the sheath region preceded by a shock, the in¬ 
tervals given for (.Pmax are not strictly speaking those of the MC 
axis. Indeed, the above values are slightly too large, for exam¬ 
ple in a well observed case, is about 10° larger for the MC 
sheath than for the FR axis ( [Janvier, Demoulin, and Dasso|20131 l. 
However, since the MC axis extension angles are not generally 
known, we report the ^^ax values estimated from CMEs in Fig¬ 
ure]^ to derive Li„t = 1.3 + 0.6 AU and Ltotai = 2.6 + 0.3 AU, 
then Ltotai ~ 2Liitt. 


5. Total Amount of Global Flux Rope Quantities 

5.1. From Local to Total Estimations 

The above results of Sections |2]|^ are applied to compute the 
magnetic axial flux F and helicity H for each FR detected at 
1 AU. The results are summarised in the distribution functions 
dLobs./dR and dHobs./dR, dependent of R. They provide the 
amount of magnetic flux and helicity per unit radius and time. 
These distributions are related to the distribution of FR number, 
dAobs./dR, as dgobs./dR = Q dNobs./dR with Q = F or H. 
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Fig. 7. Effect of the angular extension (defined in Figure[^ on (a) the curvilinear abscissa [j in AU] along the FR axis and (b) the shape of 
the FR axis for MClc set. All the curves are derived from the probability distribution [T’obsCd)] of MClc set shown in panel (c). The dashed-lines 
in panel (a) represent the case where T’obsC^) is set to be symmetric in A. This forced symmetry has a small effect on j(/l). The dotted lines in panel 
(b) represent the flux rope legs extrapolated to the Sun by a radial segment (used to computed Ltotai in Equation l|13|l). 



Fig. 8. Comparison of the results between the two analysed MC sets. The results are shown for three values of the parameter (a) The 
curvilinear abscissa [j in AU] along the FR axis (defined in Figure[TJ and (b) the shape of the deduced MC axis. The curves are derived from the 
probability distributions [T’obsfT)] shown in panel (c) with curves, rather than histogram, for comparison. The black distribution of MClc set is the 
same as in Figure |7j; and corresponds to the black, blue and green curves in panels (a) and (b). The red distribution is derived for MCuy set and 
corresponds to the red, magenta and pink curves. 


The distribution dQobs./d^ measures the local distribution of 
Q as estimated by the spacecraft over an interval of time. We are 
also interested on the total amount of these quantities crossing 
the sphere of radius D - 1 AU. In order to convert the local dis¬ 
tribution to a global distribution of FRs travelling at least up to 
1 AU, [Janvier, Demoulin, and Dasso| ( |2014a| l have estimated the 
probability to detect a FR on the sphere of radius D supposing 
a uniform distribution of FRs in longitude and in a latitude band 
+0max- This portion of sphere has a surface S^p - 47r sin Umax 
The FR extension projected on the sphere, so its apparent surface 
S, is estimated by 5 = 2RLp where Lp is the FR axial exten¬ 
sion projected radially on the sphere. The probability to detect 
this FR is Ffr = S/S^p- Then, the total distribution function is 


dQtotai/dR = dQobs./dR x 1/Pfr- This computation corrects the 
local spacecraft measurements both to estimate the total number 
of FRs launched from the Sun and to take into account the lower 
probability to detect a FR with a lower radius (as its cross section 
viewed by the spacecraft is lower). Finally, all the distributions 
are averaged over the time period of the in situ observations and 
are computed per year to be compared. 


The projected length Lp is estimated from the mean angular 
extension of CMEs close to the solar limb: i;tJinax = 30 + 15° 

1 + 0.5 AU. From the latitude 


[ Wang et al. 20111 providing Lp 
istribution of the expected solar sources, [Janvier, Demouliiu 
and Dass^ (2014a} selected flmax = 45°. The global disbibu- 


tion would be simply multiplied by a factor 0.8 (resp. 1.4) if 
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Fig. 9. Dependence of the lengths Li„,, Equation 1 12 1 , and Equation 1 13 s with the angular extension The FR has its apex located at 
1 AU and is still attached at the Sun to compute Ltotai- (a) Effect of the integration upper limit Tsup for both MCuy and MClc sets (in red and black, 
respectively), (b) The mean value and the typical range of derived from limb CMEs observed with coronagraphs (Wang et al. 2011) are used 
to estimate a mean value and a range of Lj„t and Ltotai (blue arrows) from the mean results of MCuy and MClo (blue curves). 


^max = 30° (resp. 60°) would be used instead. The results of 
the above procedure and the value ^max = 45° were back up by 
checking that the computed total number of MCs from in situ 
data matches the expected total number of MCs derived from 
coronagraph observations of CMEs (see Section [53]). 


We also study the content of helicity in small FRs in the so¬ 
lar wind. We analyse a sample of 125 small flux rope events 
presente d by [Feng, Wu, and Chao|( 2001), and also its extension 
made in |Feng et a/.| ( |2008| l (see ^Janvier, Demoulin, and Das^ 
|2014b| for details about several features of these two samples). 
The same procedure is applied to small FRs and MCs while it 
is not known whether small FRs are that broadly extended along 
their axis than MCs. Then, we simply use the same and Lp 
values for small FRs than for MCs. This choice has in fact a neg¬ 
ligible effect on the total magnetic flux and helicity estimations 
(Section [53) |, and the small FR contribution is even expected to 
be smaller as small FRs are likely to be coherent flux ropes only 
on length scales smaller than Lp used here. 


5.2. Distribution Functions 

The distribution functions can be estimated with histograms. 
However, a uniform binning of R is not suited due its large 
range of variation and the large variation of the number of FRs 
per bin. In fact the bin size should be adapted to the number 
of FRs detected in each range of R. Then, 

[and Das^ ( |2014a| l developed a technique, 
method, where the bin widths are computed to have the same 
number, A^part, of FRs in each bin in order to have a uniform sta¬ 
tistical noise across the bins. The data are first ordered by grow¬ 
ing value of R. The binning with Apart flux ropes is computed 
starting from the lower R values, then progressively shifting up¬ 
ward the bin window by one FR. This provides a smoothing of 
the fluctuations over Apart flux ropes. We use Apart = 10 as a 
compromise between decreasing the fluctuations and resolving 
the variations of the distributions. 

Figure [TO^ shows that dAtotai/d^ of MCs dominates 
dLtotai/dR of small FRs, but only by a factor 3. The MC 
contribution is peaked while the small FR contribution is almost 
independent of R (for logjyR < -2.3, the decrease is due to a 
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called the partition 


selection effect on small FR orientation, see | Janvier, Demoulin, | 
[and Dasso|2014a| l. 

Both distributions, dAtotai/dR and dHtou/dR, are maximum 
for R ss 0.13 AU with the difference that dLtotai/dR is peaked 
around this maximum while dAtotai/dR is nearly flat in the range 
0.08 <R<0.2 AU (Figure[T§i,b). For R < 0.06 AU, dAtotai/dR 
is already one order of magnitude below its plateau value so that 
the MC contribution strongly dominates for helicity. Finally, the 
slightly lower distributions for Fepping and Wu (20101, com¬ 
pared to the ones computed from the two other lists having MCs 
( |Lynch et al. 2005 [ |Feng, Wu, and Chao|2007[ ) are due to a more 
severe selection of MCs so a lower number of detected MCs by 
a factor 2. 

The distribution dAtotai/dR is a strongly decreasing function 
of R proportional to R^^ '^ in the range of small FRs (see Figure 3 
of Janvier, Demoulin, and Dasso|2014^ . A bump is present in 
this distribution for MCs, still the number of FRs is dominated 
by the small FRs (their Figure 6). By contrast, dUtotai/dJ^ and 
even more d//totai/dR distributions are larger in the MC region 
(Figure [T0^,b). The smaller FR are not enough numerous, i.e. 
dAtotai/dR is not steep enough, to balance the R^ and R^ factors 
and the less variable factors Bq and present in Equa dons ([^ 
and ([^ (Bp is on average an increasing function of R, Janvier, 
pemoulin, and Dasso|2014bl Figure 3). 

Interestingly, the distribution d//totai/dR is close to a power 
law of the radius with an exponent « 2. The dominance of 
the large scales may be a consequence of the inverse cascade 
of helicity, as found in MHD studies where helicity was found 
to be transferred from small scales to larger scales ( AlexakisT] 
|Mininni, and Pouquet|2006l and references therein). This prop¬ 
erty could he a consequence of the MHD evolution of the corona 
and even of the solar dynamo which build the coronal magnetic 
field (while flux ropes are mostly transported in the interplan¬ 
etary medium with only some erosion). The plasma composi¬ 
tion is in favor of a coronal formation ( |Feng and Wang|[2015| l, 
and other characteristics of small FRs point toward different for¬ 
mation mechanisms than MCs (see Section 5.2 in [Janvier, De-| 
jmoulin, and Dasso|2014b| l. Then, one mechanism {e.g. the tear¬ 
ing instability) cannot be put forward to explain this distribution 
of both small FRs and MCs. Since they also have similar charac¬ 
teristics at 1 AU, both the formation and the propagation process 
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Fig. 10. (a,b) Distributions of magnetic axial flux and helicity for the four studied lists of MCs and small FRs. dQ is the amount of Q(= F or H) 
in the radius interval dR per year. (c,d) Cumulative distribution functions of F and FI computed per year with an average over the time period of 
the lists (Table[^. The number of events in all curves are corrected from the apparent FR cross section projected on the sphere with 1 AU radius 
(Section[5.1|l and the summations are done from the larger to the smaller radius (Equation l|14^). FI is computed with Ltotai = 2.6 AU (Section[4.2[l. 


from an early stage in the corona to the interplanetary medium 
might play a role in transferring magnetic helicity from smaller 
scales to larger ones. We may hnd an answer to this question 
with a future study on the helicity partition at solar distances 
<1 AU. 

5 . 3 . Cumulative Functions 

We are also interested on the total amount of magnetic flux and 
helicity crossing the sphere of radius 1 AU per unit time, in order 
to have a global budget of these quantities launched by the Sun. 

We define similarly as above a global quantity by Q. We use 
below Q - 1,F,H to compute the total number, the magnetic 
flux and helicity of FRs. The amount of dQtotai in the range 
of radius dR is Q (dNtomi/dF) dR and we define the cumulative 
function of Q starting from the largest FRs since they are less 
numerous: 


Gtotal(^) 




d ltotai 

dR 


(R')dR' = 


'=X 


Q{R') ^^^^(R')dR'. 


dR 


(14) 


2totai(^) can be computed by integrating the distributions shown 
in Figure[T0^,b, or simpler by a summation of the contribution of 
each MC starting from the largest ones and applying the conver¬ 
sion factor associated to the probability to detect this FR (Pfr, 
Section mj. Indeed, the fluctuations are naturally averaged in 
a cumulative function so we show in Figure [T0);,d the results of 
this second approach which is more direct and simpler. 

The computed total number of MCs match approximately the 
number of MCs expected from the number of observed CMEs 
(| Janvier, Demoulin, and Das so |2014a|l. More precisely, the 


counts derived fromlUyncfTef a/.|(|2005|l and some CME catalogs 


are close to each other, while the same pairing is true between 
the counts derived from |Lepping and Wu| ( |2010) l and some other 
CME catalogs (see their Eigure 6). The difference of counts be¬ 
tween the pairs is about a factor 2 (Table[^ which is linked to the 
slightly different criterium used to defined both MCs and CMEs. 
This correspondence backups the above procedure which trans¬ 
forms the local measurements of Wind or ACE to global estima¬ 
tions. It shows also that the small ERs, much more numerous 
than MCs by at least a factor 10 (Table [^l, are not associated to 
CMEs. 


While the small ERs dominate by their number, they pro¬ 
vide a small contribution to Ptotai and //total (Eigur e [T0}:,d and 
Table [T]). The small ERs from the list of Feng et a/| ( |2008| l pro¬ 
vide a factor 8 lower magnetic flux and a factor 20 less mag¬ 


netic helicity per year than the MCs of the list of Lynch et al. 
( |2005] l. Both cumulative curves are indeed flat below some R 
value. The main contribution to the cumulative curves is from 
the larger MCs as expected ( Lynch et a/|2005| . Half of the con¬ 
tribution for the axial magnetic flux is from MCs with a radius 
R larger than =» 0.13 AU (with only slight variations between the 
three first lists which contain MCs, Table [^. The contribution 
for magnetic helicity is even from larger scale MCs (R > 0.14 to 
0.2 AU). 


The results from the three lists containing MCs are the clos¬ 
est for H with //total = 2.3 x 10"^^ Mx^ and at most a variation 
of 13% between lists compared to a factor 2 with the number 
of MCs. Indeed, the identification and modeling of large MCs, 
which define dominantly //total, is easier than smaller ones which 
typically have properties less contrasted, when they are com¬ 
pared to the surrounding solar-wind properties. 
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Table 1. Maximum values of the cumulative function of Ntotn, T^totai and ^totai estimated during one year. They are averages over the time period 
indicated in the first column. H is computed with Ltotai = 2.6 AU (Section [4.2[ >. The presence of small FRs is indicated by sFR and of magnetic 
clouds by MC in the third column. 


time period 

reference of ER list 

type 

^total 

number 

^total 

10^^ Mx 

^total 

10^^ Mx^ 

1995-2003 

ILynch et al. { 

2005) 

MC 

850 

31 

2600 

1995-2009 

Lepping and W 

u(|2010i 

MC 

390 

19 

2100 

1995-2001 i 

^eng, Wu, and (Jhaoj(j200 

71 MC-hsFR 

4600 

31 

2200 

1995-2005 

Feng et a/. t2008) 

sFR 

7040 

4 

130 


Table 2. FR radius [Rj, in AU] where the cumulative function of Riotai and //total (shown in Figure[To| reach half their maximum value (given 
in Table[T](. The presence of small FRs is indicated by sFR and of magnetic clouds by MC in the third column. 


time period 
1995-2003 
1995-2009 
1995-2001 
1995-2005 


reference of FR list 



Lynch et al. ( 

2005 1 

L 

epping and Wi 

I (|2010i 

Feng, Wu, and Cfi 

ao (20(1/ 


Feng et al. p 

— 


type 

/^h,N 

Rh,F 

/^h,H 

MC 

0.066 

0.13 

0.15 

MC 

0.084 

0.14 

0.20 

MCh-sER 

0.008 

0.12 

0.14 

sER 

0.005 

0.04 

0.07 


5.4. Total Magnetic Flux and Helicity from MCs and CMEs 

The amount of magnetic flux carried by MCs in one year is im¬ 
portant (Table[^. With a typical large magnetic flux of 10^^ Mx, 
but not exceptional for an active region (AR), this implies that 
all the MCs launched from the Sun carry on average the mag¬ 
netic flux of 20 to 30 ARs in one year. The amount of launched 
unsigned magnetic helicity is even more important, as follows. 
The magnetic helicity injected at the photospheric level during 
the full emerging phase of an AR is typically lO'*^ Mx^ for a flux 
of 10^^ Mx (e.g., Jeong and Chae 2007 Lim et q(.| 2007 |Tian 


and Alexander |2068[ Jing et al. 2012| l. This implies that MCs 

carry per year the magnetic helicity of about 200 emerging ARs. 

The above estimations of Ftotai and //total are computed only 
for the estimated total number of MCs crossing 1 AU in one year. 
However, MCs are detected in about only one third of ICMEs on 
average over the solar cycle (e.g., Richardson and Cane||2010| 
and references therein). In about one third of ICMEs, a mag¬ 
netic held rotation is detected but it is not coherent enough or 
a low enough proton temperature is not found, so they are not 
classified as MCs but as cloud-like events ( |Lepping, Wu, and] 
|Berdichevsky|2005| . If the lower detection rate of MCs is sim¬ 
ply due to the spacecraft passing on the side or missing the flux 
rope, as advocated by |Jian et a/.| ( |200^ , so if all ICMEs have a 
flux rope inside (as recent results point out: |Gopalswamy et aL\ 
|2013[|Makela et a/|2013[ ), the amounts of /^totai and //total shown 
in Eigure 10 and Table [Ijneed to be multiplied by a factor 3. 

On the contrary, the amount of estimated helicity could be 
lower because some MCs are likely no longer attached to the 
Sun when observed at 1 AU. Not taking into account the length 
of MC legs decreases the length and the helicity estimated by a 
factor 2 (Section [4.2| l. Next, the erosion of the ER by reconnec¬ 
tion with the solar-wind magnetic held is not taken into account 
in the above lists and deduced results. The analysed MCs contain 
the ER remaining intact at 1 AU, and a part of the reconnected 
flux (called a back region by |Dasso et al.\2QQ6\ . Then, the above 
flux and helicity estimates are in between the ones from the re¬ 
maining ERs at 1 AU and the ERs launched from the Sun. Since 
the average amount of reconnected flux is large, about 40 % of 
the total azimuthal magnetic flux, the above helicity estimate is 
expected to be a factor around 2 too low for the helicity launched 
by the Sun. Einally, the ER cross section is typically flat by a 
factor 2 to 3 on average ( Demoulin, Dasso, and Janvier||2013l l. 
Compared to the cylindrical model used to fit the data, this in¬ 


troduces an underestimation of the helicity by a factor slightly 
below 2 to 3 (helicity is near proportional to the aspect ratio for 
bja > 2, see Eigure 8a of |Demoulin and Dasso|2009b) l. To sum¬ 
marise, the over estimation of helicity, implied by supposing all 
MCs observed at 1 AU as attached to the Sun, is likely to be 
over-compensated by the other factors described above. Then, 
our estimates of magnetic helicity (Eigure [TO] Tables flp] ! are ex¬ 
pected to be underestimated by at least a factor 2, and plausibly 
a factor 6 if non-MC CMEs carry the same amount of helicity 
than MCs. 


5.5. Helicity Budget over a Solar Cycle 

Below we estimate the total unsigned helicity that leaves the Sun 
over a solar cycle. The list of |Lepping and Wu| ( |2010| l covers 
most of a solar cycle 23. While the two other lists are more re¬ 
stricted in time. Table [T] they provide comparable values of he¬ 
licity transported per year. Then, we convert these results to the 
total amount of unsigned magnetic helicity transported by MCs 
during one solar cycle by both hemisphere (assuming a constant 
mean helicity per year): //mc, cycle ~ 2.5 x 10"^® Mx^. This is 


one order of magnitude larger than the estimation of Bieber and 
Rust ( 1995|l from solar cycles 20-22, and a factor 2.5 larger than 


the one of DeVore| ( [2000| l from solar cycle 21. However, in this 
last case the difference is mostly due to difference of axis length 
used: 0.5 AU for DeVore ( |2000| l compare to 2.6 AU here (Sec- 
tion |4.2| i. Then, within a factor 2, our results agree with the order 
of magnitude estimated by DeVore ( 2000) 1. 

Before comparing our result to the solar source estimate, it 
is worth to compare it to other solar phenomena to appreciate 
its magnitude, //mc, cycle is three orders of magnitude larger than 
the total helicity injected in the quiet Sun ( |Welsch and Longcope) 
2003|l and one order of magnitude lower than the total helicity 


injected in open held of coronal holes ( Berger and Ruzmaikin) 
|2000) l while the total unsigned magnetic flux involved have com¬ 
parable magnitudes. Next, the analytical expression for the mag¬ 
netic helicity contained in a simplified Parkerian solar wind for a 
period of a solar rotation was computed by |Bieber, Evenson, and) 
Matthaeus ( |1987| l (see their Equation (8)). Erom this expression, 
an helicity near 7x10"^^ Mx^ is obtained for a complete solar cy 


cle. It is a value co mparable to the solar estimation of Berger and 
( 2000|l with a constant open flux ^ 4x 10^^ Mx per 


Ruzmaikin 


magnetic polarities. Then, the ejection of ERs from the Sun is an 
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efficient mechanism to eject magnetic helicity, but less efficient 
than the direct solar rotation twisting the open flux. 

Next, the MCs are not related to the quiet Sun nor to the solar 
open flux, but to the solar dynamo and ARs, as follows. Differen¬ 
tial rotation in the convection zone creates an opposite amount 
of magnetic helicity in each solar hemisphere. The amount of 
unsigned magnetic helicity create d during a solar cycle is about 
ss 5 X 10"^^ Mx^ for solar cycle 22 (Berger and Ruzmaikin 2000 1 
so only a factor 2 larger than the above estimate for MCs. The 
amount of helicity created by the a effect is more difficult to es¬ 
timate, but they argued that the amount is comparable or larger 
than the amount provided by differential rotation. Then the solar 
dynamo is able to create 4 or more times the amount of unsigned 
helicity found in MCs. If most CMEs have a FR, and taking 
into account the flatness of the FR cross-section, the amount of 
unsigned helicity transported by CMFs is similar to the one pro¬ 
duced by the solar dynamo. 

The magnetic held, amplified by the global dynamo, emerges 
mostly in ARs. Improved local correlation tracking methods 
have been developed to derive the photospheric velocities. From 
these measurements magnetic helicity fluxes are derived (e.g. see 
the review of Demoulin and Pariat|20091 l. The largest input of 
helicity in the solar atmosphere is detected during the emergence 
of ARs. //MC.cyde is about a factor 4 larger than the amount of un¬ 
signed helicity injection, ^ 0.6 x 10^^ Mx^, found by Georgoulis 


\et al.\ ( |2009| l in emerging ARs over solar cycle 23. However, 
//mc, cycle IS about the value found by Yang and Zhang (2012|l, 


3.3 X 10"^^ Mx^, and half the value found by Zhang and Yang 
5 X 10^® Mx^ for AR emergence during solar cycle 23. 


( |2013| l 


We conclude that, within the present uncertainties of mag¬ 
netic helicity estimations. Section [5^ the amount of magnetic 
helicity sent in MCs/CMFs is compatible both with the amounts 
of helicity built by the solar dynamo and measured in emerging 
ARs. 


6. Conclusion 


The in situ measurements only provide local information of the 
physical parameters along the spacecraft trajectory. These data 
are typically fitted by a FR model to estimate the FR properties 
in a local 2D cut orthogonal to the FR axis. Since explorations of 
the same MC by several spacecraft are rare, we used a statistical 
approach to derive the generic properties of MCs. Indeed, dif¬ 
ferent MCs are crossed at different location along their axis, pro¬ 
viding statistical information along the axis. Janvier, Demoulin, 
[and Dass^ ( |2013| ) developed this new technique and they found 
the generic shape of MC axis only parametrised by the angular 
extension ^max- The location of the spacecraft crossing along the 
axis is related by the location angle A which is available for each 
MC from the axis direction determined by fitting a FR model to 
the in situ data. 

In the present study, we further develop this statistical 
method by first analysing how FR quantities vary along the FR 
axis and second by computing the axis length. We And no de¬ 
pendence of the FR radius, held strength and twist along the FR 
axis in a broad range around its apex, within |/l| < 50°, for MCs 
observed at 1 AU. The variation found for larger |/l| values are 
interpreted as a bias introduced when the spacecraft cross the FR 
legs and explore a significant portion of the FR along its curved 
axis while the fitted model has a straight axis. We And a mean 
axis shape nearly symmetrical on both sides of the apex, with a 
very little asymmetry qualitatively consistent with the expected 
deformation given by the solar rotation. Next, we derive the 


length of the generic axis: 1.3 + 0.6 AU where the uncertainty is 
derived from the range of ^max observed for limb CMFs (Wang 
\et a/.|20lTl l. If the FR is still attached to the Sun, the minimum 
length is 2.6 + 0.3 AU. The results allow the transformation of 
the local estimation of magnetic helicity, per unit length, to the 
total helicity of the FR. 


The above results were applied to four lists of events: two 
with only MCs, one with MCs and small FRs and one with only 
small FRs. While the small FRs largely dominate MCs in num¬ 
ber (taking into account the probability to detect a FR), they have 
a much lower contribution (at least by a factor 10) than MCs for 
magnetic flux and helicity. Indeed, MCs transport an important 
amount of magnetic flux and helicity when estimated over the 
full 1 AU sphere, as follows. During one year, MCs carry a 
magnetic flux Fyic ~ 27 x 10^^ Mx^ and an unsigned magnetic 
helicity Hmc ~ 2.3 x 10^^ Mx^. These are equivalent of the flux 
contained in about 25 large ARs and the equivalent of helicity 
injected in 200 emerging and large ARs (with a magnetic flux 
of ss 10^^ Mx). If all ICMFs possess a FR component, these 
numbers have to be multiplied by about a factor 3. Finally, the 
amount of unsigned magnetic helicity carried away from the Sun 
by MCs during a solar cycle is comparable to the amount esti¬ 
mated for the solar dynamo and to the one measured in emerging 
ARs. 


While we have improved the helicity estimation in MCs by 
analysing the FR parameter dependence along the axis and es¬ 
timating the FR length, there are still a number of issues which 
require improvements. A first one is the local FR model used to 
fit the magnetic data. Does this model characterise well enough 
MCs? Comparable helicity values were found with different 
models ( |Gulisano et a/.|[2005| l, still it would be worth to do 
broader explorations both in terms of MCs and models, espe¬ 
cially since doubts on the relevance of the Fundquist’s model 


have recently come out ( 

Hu et al. 

2014)1. Second, FRs erode as 

they propagate in the so 

ar wind ( 

Dasso et al. ||2006[ [Ruffenach 


\et a/.|[20T5] l. A deeper analysis of present data would allow to 
estimate both the helicity remaining in the FR at 1 AU and the 
one present before reconnection. Third, many FRs do not have 
a circular cross section, so that an effort to fit e.g. an elliptical 
model to the magnetic data would improve the helicity estima¬ 
tion. Finally, the solar helicity budget can be determined over 
the same time interval during which the MC helicity budget is 
studied. We conclude that there is a real potential to further im¬ 
prove our knowledge of MCs and in particular the solar magnetic 
helicity budget. 


Appendix A: Mean Twist of a Flux Rope 


We consider in this section any magnetic held with a cylindri¬ 
cal symmetry, so B(r) = B^ir) Ca + B^{r) Cz, where Ba, B^ are 
the azimuthal and axial components depending only on the ra¬ 
dial coordinate r. For this FR configuration, the magnetic he¬ 
licity is linked to the amount of turns [n(r)] per unit length as 
shown below. This is a concrete application of the more gen¬ 
eral expression of Equation (12) of Berger and Field] ( 1984i with 
poloidal/toroidal decomposition of a magnetic held. However, 
the derivation rather follows the work of [Berger] ( |2003| l with the 
mutual helicity of “open” held (here B^) and “closed” held (here 
Ba) which is gauge invariant. 

The magnetic helicity. Equation (j^, involves the vector po¬ 
tential component Aa. Since V x A = B, Aa is linked to Bz as 
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(Equation (2) of |Dasso et a(.|2005b| l 
r , , , F{r) 

rAa(r)= r'B,(r')Ar'^ (A.l) 

Jo 27r 


where F(r) is the axial magnetic flux within the circle of radius 
r. Writing the held line equations, the number of turns n{r) is a 
function B components as 


n{r) = 


BJr) 

In r BJr) ' 


(A.2) 


Inserting Equation ( |A.l| l in the left equality of Equation Q and 
replacing B^ir) using Equation (|A.2|l implies 


H^2L f 
Jo 


F(r) n(r) 2n r B^{r)Ar 


=2L r 

Jo 


F(r) n{r) 


AF{r) 

dr 


dr, 


(A.3) 


where L is the length and R the radius of the ER. Then, FI is 
rewritten as an integral on the axial flux F as 

ff^2L f n(F')F'dF'. (A.4) 

Jo 

When n is independent of the radius r, so of F(r), Equation ( |A.4| l 
reduces to 

H^nF^L. (A.5) 

More generally, one defines with Equation 0 as a flux 
weighted mean of the number of turns per unit length. Then, the 
magnetic helicity of a ER is always of the form of Equation ( |A.5| l 
with n replaced by rit, as written in Equation Q. 
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